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PREFACE 


This  report  is  part  of  RAND's  continuing  theoretical  study  of 
the  properties  of  high  temperature  plasmas.  The  results  are  of 
interest  for  application  to  the  determination  of  densities  and 
temperatures  of  high  altitude  plasmas. 

Several  installations  are  under  construction  to  measure  the 
backscatter  from  powerful  radar  beams  directed  vertically  upward, 
and  it  is  hoped  that  incoherent  scattering  of  electromagnetic  waves 
from  a  laboratory  plasma  may  be  observed. 

These  experiments  and  associated  theoretical  studies  will 
lead  to  a  better  understanding  of  geophysics  and  electromagnetic 
wave  sciences. 


SUMMARY 


The  incoherent  scattering  of  radio  waves  from  a  hot  plasma  is 
confuted  using  the  diagrammatic  techniques  of  quantum  electrodynamics 
This  method  simplifies  previously  obtained  results  and  permits  ex¬ 


tensions  to  include  the  effects  of  close  collisions* 
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I.  IWTRODUCTIOW 

Several  years  ago  W.  E.  Gordon^ suggested  that  the  weak  hut 
measurable  incoherent  scattering  of  radio  waves  from  electrons  at 
high  altitudes  would  provide  information  about  their  density  and 
temperature  out  to  a  distance  of  several  thousand  kilometers*  At 

2  1/2 

frequencies  well  above  the  electron  plasma  frequency  ®  V®)  ^ 

the  ionosphere  is  essentially  transparent.  Radio  waves  are  then 
scattered  by  charge  fluctuations^  and  the  scattered  power  is  proportional 
to  the  number  of  particles.  If  the  particles  did  not  interact  eimong 
themselves  one  would  obtain  for  the  scattering  cross-section  N  times 
the  familiar  individual  Doppler-spread  cross-section: 

-  5  <i)' 

£  (»,k)  •  »  r/  e  (1.1) 

2  2 

where  N  is  the  number  of  scatterers,  r  =  e  /m  c  is  the  classical 

o 

electron  radius,  cj  »  ((1)5  -  •‘laVfip  is  the  frequency  shift,  and 
k  =  (k.jj  -  Is  the  (vector)  change  in  wave  nuniber  In  units  of 

2  1/2 

the  Debye  wave  number  ^  ®  n/kT)  •  The  cross-section  for 

2 

backscatter  is  then  or,  ^  ^  =  H  r 

total  o 

(2) 

An  experiment  was  performed  by  K«  W.  Bowles*  He  observed 
the  incoherent  scattering  but  found  that  the  frequency  spread  of 
the  retixmed  signal  was  much  narrower  than  exi)ected,  and  he  proposed^ 
that  the  spread  of  the  returned  signal  was  characteristic  of  the  ion 

velocity  and  not  the  electron  velocity.  Subsequent  theoretical  in- 

fU)  (5)  (6-9) 

vestigations  by  Salpeter,'  '  Dougherty  and  Parley,'^'-  and  others' 

have  confirmed  this  conjecture  and  have  presented  a  more  detailed 
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picture  of  the  scattered  radiation.  The  shape  of  the  observed  signal 
is  determined  largely  by  collective  effects. 

We  shall  show  here  how  the  standard  results  are  obtained  in  a 
simple  and  transparent  way  from  a  diagrammatic  theory  of  many-body 
electromagnetic  interact  ions.  This  point  of  view  enables  us  to 

include  the  effects  of  (close)  collisions,  which  have  been  neglected 
by  previous  authors. 
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II.  FORMUIATKW 


The  description  of  Individual  electromagnetic  processes  In  terms 
of  Feynman  diagrams  Is  at  once  the  clearest,  the  most  elegant,  and 
the  most  convenient  for  the  purpose  of  calculation.  We  shall  utilize 
this  method  for  the  problem  at  hand:  the  Incoherent  scattering  of 
photons  from  a  weakly-coupled,  hlgh-tempe nature  electron-ion  plasma 
in  thermodynamic  equilibrium. 

The  total  scattering  rate  Is  given  by  the  familiar  Golden  Rule 


Ftotal 


2 


<  2; 


(2.1) 


(2ii)V(P2-Pj^+  *  *  \)  2*6(E2-E^+  *  -fi  . 


The  subscripts  1  and  2  refer  to  the  initial  and  final  states  of  the 

particles.  The  subscripts  a  and  b  refer  to  the  initial  and  final 

states  of  the  radiation  (k  is  the  initial  wave  number,  e  is  the 

a  a 

initial  polarization).  In  Eq.  (2.1)  we  sum  over  all  final  states  and 
average  over  all  initial  states.  The  correct  average  over  the  initial 
particles  states  is  obtained  by  using  a  suitable  Gibbs  factor, 

e"  (2.2) 

where  P  ^  kT,  u  Is  the  chemical  potential,  Nj^ls  the  number  of  parti¬ 
cles  In  the  state  1,  and  exp  pO  Is  the  normalizing  factor  fixed  by 


Trpj^  «  1 
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We  are,  however,  Interested  not  in  the  total  transition  rate  but 
in  the  rate  to  all  particle  states  with  the  final  photon  state  fixed. 
We  denote  this  by  f‘(k,u>).  From  Eq.  (2.1)  one  has 


r(k,u))  -  Im|  1;  \  > 

e.  e  2  1  ® 

d'  a 


(2.3) 


(2*)^  6^(P2-P^+  tk)  2x6(E2-  E^^  -  ^) 


To  obtain  the  amplitude  M  we  simply  draw  all  possible  modes 
for  scattering  a  photon  in  state  (k^,e^)  to  a  state  (k^,e^).  We 
can  make  a  number  of  good  approximations.  In  the  non-relativistic 
limit  the  diagrams  with  single  photon  vertices  can  be  neglected  in 
comparison  with  diagrams  with  a  double  photon  vex*tex.  Furthermore^ 
we  shall  neglect  terms  of  order  a  «  m/u,  the  ratio  of  the  electron 
and  ion  masses,  so  it  is  not  necessary  to  consider  the  interaction 
of  the  radiation  with  the  ions.  In  Fig.  la  we  show  the  slnq>lest 
interaction  of  Interest  here.  Other  possible  processes  are  shown 
in  Figs.  1-3*  The  braided  line  represents  a  screened  Coulomb  inter¬ 
action. 

The  amplitude  for  any  process  is  readily  computed  from  the 

corresponding  diagram  by  use  of  the  calculating  rules  presented  in 
•* 

Ref.  10.  For  example,  the  amplitude  for  the  double  photon-electron 

veirtex  shown  in  Fig.  la  is  given  in  Ref.  10  to  be  (4ite  /m)^.%^. 

The  rate  for  this  process  is  then 
* 

The  amplitude  may  be  infinite  when  the  high  frequency  expansion 
for  that  process  is  not  valid.  It  is  then  necessary  to  use  a  modified 
perturbation  expansion.  That  is,  one  ssist  first  siss  certsdn  classes 
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r(k,u.)-i^  j 


_s,r  1 1 

(ac)’  (a.)’  "  8?“* 


Ulte  ^ 

■■■  e.  .e 
m  b  a 


(2.4) 


(2it)V(p2-p^  +  &  k)  2jc  +  &  «) 
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We  shall  deal  only  with  unpolarized  radiatioi^  so  we  siun  and 
average  over  the  polarizations. 

k  I  -5(1+  cos^O)  (2.5) 

®b'®a 

where  backseat tcred  radiation  one  readily 

obtains  the  result  given  in  Section  1  for  non*lnteracting  particles. 

In  the  following  section  we  shall  discuss  in  detail  the  many- 
body  processes  that  produce  the  central  characteristic  resonance  in 
the  scattered  radiation. 


of  diagrams  to  remove  the  divergence  at  low  frequency. 

Notice  also  that  the  rules  for  the  distribution  functions  are 
somewhat  different  (and  simpler)  her^  since  %ie  want  only  the  forward 
process • 
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III>  STAI^ARD  RESULTS  BY  THE  DIAGRAM  METHOD 


The  scattering  due  to  the  process  shown  in  Fig.  la  produces  a 
rather  broad  line.  This  is  a  property  of  the  one-electron  pair  final 
state  and  so  is  also  true  of  the  process  shown  in  Fig.  lb.  Since  we 
are  interested  here  mainly  in  the  sharp  central  line  we  shall  disre¬ 
gard  these  processes  and  concentrate  on  the  process  shown  in  Fig.  Ic 
whose  ion-pair  final  state  produces  a  narrow  line. 

Direct  scattering  of  radiation  by  an  ion,  analogous  to  the 

electron  process  in  Fig.  la,  also  produces  a  narrow  line,  but  is 

less  important  than  the  process  shown  in  Fig.  Ic.  There  are  two 

reasons  for  this:  (l)  the  direct  photon-ion  scattering  has  an  extra 

4  2 

factor  of  a  =  (hi/m)  ,  and  (2)  the  matrix  element  for  the  indirect 
process  (Fig.  Ic)  has  a  resonance  at  the  low  frequency  acoustic  mode. 

The  physical  interpretation  of  Fig.  Ic  is  that  the  photon  (wavy 
line)  is  scattered  by  an  electron  (the  bubble)  in  the  screening  cloud 
which  accompanies  the  ion  (heavy  line)  and  interacts  with  it  by  means 
of  a  screened  interaction  (braided  line).  The  amplitude  for  this 
process  is  given  by  (amplitude  for  double  photon-electron  vertex)  x 
(amplitude  for  electron  bubble)  x  (amplitude  for  screened  Coulomb 
interaction),  or  in  the  notation  of  Ref.  10, 


®b'V  (-  i  — 


k  +  Q(k,a)) 


■) 

/  o 


2  Q. 


-  ®b-®a(5-l) 

k  +  Q 


where  \ 


The  function  Q(k,a))  represents  the  entire  proper 


I>olarlzatlon  part.  For  the  purposes  of  this  section  we  want  Q  only 
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to  lowest  order,  Q  ^  >  the  sum  of  the  electron  and  Ion 

hubbies. 

The  scattering  rate  at  0  *  n  is  then  given  by 


r(k,„).r 

'•  J  8,=  “a 


^3/2 


1  2 


‘^e 

?  k^+Q 


(2«)5b^(P2-Pj^  +  fi  k)  2*  6(i  a\‘^  -  i  a\‘'  +  ft  u) 

t ^  (2 


,12  2  12  2 


n  r 


Ik^  +  q!^  ViJt  ock 


1  i/ 

1  ■"  ?? 
e 


(3.2) 


If  we  take  the  values  of  and  Q  In  the  random  phase  approximation 
then  It  Is  easy  to  show  that  this  Is  Just  the  result  obtained  by 
previous  authors  (by  rather  less  direct  methods).  The  subscript  aonfj^ 
refers  to  the  acoustic  mode.  The  analogous  result  can  be  obtained 
near  the  plasma  resonance  by  putting  a  >  1  in  Eq.  (3*2).  It  should 
be  noted  that  these  results  are  accurate  only  near  the  respective 
resonances.  It  is  easy  to  write  down  a  result  vdilch  is  valid  over 
the  whole  range  of  frequency  but  it  is  more  conqpllcated.  In  Fig.  4 
we  show  several  plots  of  the  scattering  cross-section  in  the  random 
phase  approximation  for  several  values  of  k  .  From  Eqs.  (6.12)  and 
(6.13)  of  Ref*  10  we  have  (with  Z  «  u/ok) 


qJ(z) 


_  1  2  2  12 

1-Ze^  Jdte^  +i  (*/2)^/^  Z  e  ^ 


(3.5) 


and  -  ^^(az) 
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In  general,  then,  we  can  describe  the  result  as  follows.  When 
collective  effects  are  unimportant,  that  is,  when  the  wave  nvunber 
k  =  is  large,  the  scattering  is  mainly  through  the  mode  shown 

in  Fig.  la.  As  the  wave  number  transfer  k  decreases  collective  effects 
become  increasingly  important  and  processes  of  the  type  shown  in 
Figs,  lb  and  Ic  dominate  the  scattering.  The  denominator  in  Eq.  (5*2) 
is  the  absolute  square  of  the  dielectric  function  of  the  medium. 

The  scattering  is  thus  enhanced  at  the  resonant  frequences  of  the 
medium.  There  are  two  such  frequences:  the  high  frequency  plasma 
mode  and  the  low  frequency  acoustic  mode.  In  the  center  of  the 
scattered  line,  where  the  scattering  is  sharply  enhanced  by 

the  acoustic  resonances  on  each  side  of  cOg.  Since  this  mode  is  fairly 
strongly  damped,  at  least  when  the  electron  and  ion  temperatures  are 
equal,  the  two  resonances  run  together  and  the  peaks  are  Just  barely 
distinguishable.  At  about  ~  i  sharp  resonances 

corresponding  to  the  electron  plasma  frequency. 
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IV.  THE  EFFECT  OF  cnT.T.TRTrwR 


To  obtain  the  lowest  order  effect  of  (close)  collisions  ve  must 
compute  r (k,u)  beyond  the  random  phase  approximation.  In  order  to 
correct  the  shape  of  the  entire  spectrum  of  the  scattered  radiation 
It  Is  necessary  to  have  expressions  for  the  collision  corrections 
(or,  the  conductivity)  \dilch  are  valid  for  the  entire  range  of  fre¬ 
quencies,  Including  very  low  frequencies.  Such  expressions  are  not 

* 

available  in  the  literature  and  the  only  results  which  have 
appeared^ are  valid  only  at  high  frequencies,  that  is, 
when  the  frequency  taken  up  by  the  system  Is  greater  than  some 
appropriate  collision  frequency.  Nevertheless,  one  finds  that  for 
certain  cases  the  high-frequency  formulas  apply  for  a  part  of  the 
region  of  Interest.  They  always  apply  sufficiently  far  from  the 
center  of  the  line. 

We  are  Interested  mainly  In  those  processes  which  broaden  the 
final  states  in  Fig.  Ic  so  that  the  infinitely  sharp  energy  conserving 
6-function  in  Eq.  (3*2)  becosies  a  finite  resonance  function 


1 

it 


ii  r 


(*U)  + 


2  2 

gapg 


■  2  “ 


‘^1 


(4.1) 


The  width  is  the  transpoirt  collision  frequency  for  a  particle 
of  momentum  p.  It  is  the  svim  of  all  partial  widths  for  the  various 
collision  processes  which  alter  the  ion  motion. 


One  of  the  authors  (D.  F.  DuBols)  has  recently  obtained  such  a 
result  for  the  conductivity  of  an  electron-ion  plasma.  (To  be 
published.) 
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We  flnd^  therefore^  that  the  Important  processes  are  the  ones 

shown  in  Figs#  2a  and  2b#  The  high-frequency  rates  are  readily  com- 

* 

puted  from  the  results  of  Ref#  10, 

To  obtain  a  generally  valid  result  one  must  in  effect  consider 
all  iterations  of  the  scattering  processes  in  Fig#  2.  At  low 
frequencies  the  result  of  this  iteration  is  a  Boltzmann- like  equation 

whose  solution  yields  ^p# 

Other  collision  processes  are  shown  in  Fig#  3#  The  process 
shown  in  Fig#  3a  is  tinimportant  because  it  lacks  an  ion  resonance# 
Note,  however,  that  both  processes  are  of  the  same  order  in  the 
interaction  parameter  ^  because  there  is  a  factor  ^  ^  for  each  closed 
loop#  The  process  shown  in  Fig#  3b  vanishes  in  the  limit  of  small 
frequenc ies# 


it 

Strictly  speaking,  the  following  results  are  valid  only  at  very 

high  temperatures,  say  above  10^  °K,  when  the  Born  approximation  is 
valid.  At  lower  temperatures  one  expects  a  change  only  in  the  argu¬ 
ment  of  the  logarithm  corresponding  to  the  change  in  the  distance  of 
minimum  approach  from  the  thermal  deBroglie  wave  length  to  e^/kl# 
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V>  DISCUSSION 


A  number  of  problems  remain  to  be  solved  in  order  that  collision 
corrections  be  available  for  the  scattering  cross-section  for  a  wide 
range  of  parameters.  The  main  problem  is  to  obtain  a  formula  for 
the  conductivity  of  a  plasma  which  is  correct  for  all  frequencies. 

If  this  were  available  it  would  be  possible  to  extend  our  results 
to  the  center  of  the  line. 

For  application  to  scattering  from  the  ionosphere  it  is  important 

to  include  the  effect  of  a  magnetic  field.  This  depends,  of  coxirse, 

on  the  frequency  of  the  radiation  .  The  neglect  of  the  magnetic 

,  (^) 

field  is  not  a  bad  approximation  at,  say,  400  Me. 

There  are  two  interesting  extensions  of  this  problem.  The  first 

is  to  permit  different  tenqperatures  for  the  electrons  and  the 
(4  9) 

ions.'  When  the  electrons  are  much  hotter  than  the  ions  the 

acoustic  mode  is  very  weakly  damped  and  the  two  peaks  in  the  central 

r9) 

line  become  quite  sharp.  The  other  extension'  is  to  allow  a 
relative  velocity  between  the  two  components.  This  can  further 
sharpen  the  resonance,  with  a  very  large  increase  in  cross-section, 
until  at  some  critical  velocity  the  entire  system  becomes  unstable. 

It  would  be  Interesting  to  find  the  effect  of  collisions  on  these 


instabilities 
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Appendix 

CONNECTION  WITH  FAMILIAR  METHODS 
In  this  appendix  we  present  the  problem  in  an  alternate  formula¬ 
tion  which  makes  closer  contact  with  more  familiar  methods.  We  obtain 
the  well-known  result  that  the  scattering  cross-section  is  propor¬ 
tional  to  the  imaginary  part  of  a  density  correlation  function. 

The  following  account  is  self-contained^  but  brief.  The  in¬ 
terested  reader  will  find  a  more  detailed  discussion  of  electro¬ 
magnetic  interactions  in  plasmas  in  Ref.  10. 

The  part  of  the  interaction  Hamiltonian  which  is  significant 
for  this  problem  is 

H,  (t)  ■ i  d^x  n(x,t)  A(x,t)  •  A(x,t)  .  (A.l) 

^  2mc 

The  electromagnetic  field  is  described  by  the  (time-dependent) 
operator  A(x,t)  and  n(x,t)  is  the  electron  density  operator. 

The  amplitude  for  a  transision  1  2  is  then  proportional  to 

J  dt  <  2;  1  H^(t)  1  1;  k^,e^  > 

2 

■  — T7?  ^  ^  ^  I  n(o,o)|  1  )  (A. 2) 

'  b  a 

(2n)3‘^f2  -  '"S  -  -  El  +  h  “-i,  -  h  V- 

One  then  imnedlately  obtains 
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r(k,»)  Pli;  r  (%■%>' 

2  1  ^ 


<2|n(o,o)ll> 


(2jt)V(p2-P3^ 


(A.5) 


ConQ)arison  of  Eq.  (A. 5)  with  Eq.  (4.12)  of  Ref.  10  shows  that  the 

transition  rate  is  proportional  to  the  imaginary  part  of  a  density 

•)(.  ^ 

correlation  function  r(k  ,0)).  We  can  rewrite  Eq.  (A. 3)  in  the  form: 


r- /,  \  ^o  ^  .  _2  Im  7T  (k,(o) 

I  (k.w)  =  — —  (1  +  cos  0)  1  ‘  • 

2j(u)' 

a  1-e 


(A-4) 


At  high  temperatures  we  can  put  l-e^^saha). 

The  reader  is  urged  to  consult  Ref.  10  for  a  precise  definition 
ofX]'^(k,cj)  and  the  details  of  the  diagrajnmatic  analysis  of  this 
function.  The  open-diagram  techniques  lead  directly  to  the  results 
of  Section  III. 

In  Fig.  5  show  the  diagrammatic  description  of  TT,  the  proper 
polarization  part  Q,  and  the  screened  interaction  V*  The  correlation 
function  TT  can  be  expressed  in  terras  of 


4ir  7r(k,u))  «  Q  -  Q  -■  p  Q  , 

1+k*  Q 


(A.5) 


and#  if  we  drop  terras  which  vanish  with  a  =  m/M  and  put  in  for  Q  the 

expression  obtained  in  the  random  phase  approximatiot^  we  immediately 

get  the  well-known  result  for  p (k#u)  vdilch  has  been  obtained  by 

(4-9) 

previous  authors . ' 


*The  function  differs  from  the  TT^j  of  Ref.  10  only  in  that  the 
current  densities  Jjl(o)  and  Jj(o)  are  replaced  by  charge  densities. 


Fig.  1b 
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Fig.  1  c 
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Photon  scattering  by  a  free  electron  (la),  by  an 
electron  In  the  cloud  of  another  electron  (lb), 
and  by  an  electron  In  the  cloud  of  an  Ion  (lc)> 


Resonance  scattering  to  a  two  ion-pair  final 
state  (2a)  and  to  a  one  ion-pair,  one  electron- 
pair  final  state  (2b)» 


.16- 


Fig.  3a 


Other  collision  processes 
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Fig.  4 


Scattering  line  shape  in  the  random  ptose 
approximation  for  several  values  of  k^. 


Fig.  5 


Dlagraimatlc  description  of  ir»  Q>  and  V 
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